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The damping mechanisms of the A (1232) resonance in nuclei are studied by analyzing the 
quasi-free decay reactions ^^C(7r"'',7r^p)^^B and ^^C("^He,t7r^p)^^B and the 2p emission reactions 
^^C(7r+,pp)^°B and ^^C(^He,tpp)^°B. The coincidence cross sections are calculated within the frame- 
work of the isobar-hole model. It is found that the 2p emission process induced by the decay of the 
A resonance in the nucleus can be consistently described by a tt + p -f p' model for the A-|-N— > N-l-N 
. . . decay interaction. 

o 

(D . I. INTRODUCTION 

Q 

J — ■ Recently, new information on the A(1232)-propagation in the nucleus has been obtained from a coherent pion decay 
^—1 \ experiment [|l|,|j and from a (p, n) spin-flip transfer experiment . In the flrst experiment the ^^C('^IIe,t7r+)^^C(g.s.) 
reaction was used to measure the isovector spin- longitudinal [S ■ qT) response function in the A resonance region. In 
the second experiment the spin observables of the (p, n) reaction were used to decompose the charge exchange cross 
section in the A resonance region into its spin-longitudinal (LO) and spin-transverse (TR) components. Similar to the 
TT-nucleus total cross section data m|,0,^ the LO cross sections of both reactions show a substantial downward energy 
, shift of the A resonance in nuclei, as compared to the proton target. From a consistent A-hole model analysis of pion 
p\| ' and photon scattering, and charge exchange reactions [^|j|,|l^,|ll|] it is found that a large part of the observed shift is 
T-H , due to a nuclear medium effect on the LO response function. The medium effect is caused by the strongly attractive, 
energy dependent A-particle - nucleon-hole residual interaction VAAf.AAr. In Refs. j|J|,|l^,^,^ it was shown that 
\ Van, AN can be well described by the n + p + g' model ( and references therein) . The strong attraction of the 
(~| . TT-exchange in the LO channel produces a collective pion mode at excitation energies of ~ 250 MeV in the laboratory 
' frame. The collectivity shifts the LO response function down in energy by 60 MeV relative to the spin-transverse 
' {S X qT) response function. Other, smaller effects come from A conversion processes, such as AA^ N + N [p|,p^, 
^ ■ and from projectile excitation p^ , p^ . 

P5 , In the present paper we apply the A-hole model of Refs. to the calculation of the damping of the collective 

^ ' pion mode in the nucleus. The major decay channels are the coherent pion decay, the quasi-free decay, and the 
, 2p emission. While various calculations for the coherent pion decay were published already in Refs. |Jll|,|l|,|l|,0 
■ we give here the results for the quasi- free (7r"'"p) decay and the 2p emission. Since the coupling interaction for the 
H ' quasi-free decay is known we can use this process to study the distortion effects on the outgoing pion and proton 
. . 1 wave functions. For the 2p emission process we assume a. n + p + g' interaction. We show that the 2p emission in the 
A resonance energy region is dominated by the zero-range Landau-Migdal interaction, the strength of which can be 
exctracted from the data. Both the pion induced reactions and the charge exchange reactions are well reproduced by 
calculations with a Landau-Migdal parameter in the range of « 0.25 — 0.35. 



II. THEORY 



In this section we describe the A-hole model used in the analysis of the experimental data. The formalism and 
the methods of calculation were presented already in recent papers [^|JTl|] . In the present paper we discuss only those 
formulas which are connected with the quasi-free decay and the 2p emission. 
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A. The excitation processes 



We start our formulation by writing down the fivefold differential cross section for the charge exchange reaction 
A + a^{B + A)+b^C + c + d + b'mthc LAB frame {Ea = m^, pA = 0) 

d^cr 11 TUa rrib nic rud mc Pb Ef, p^ E^ Pd E^ Ec 



dEbdQ.iydEcdQ.cdVLd {2ttY \ v^ei \ Ea Eb E^ Ed Ec \ (mA + i^c.e. - E^) - -^(pa cos 9ad - Pb cos dbd - Pc cos 0cd) \ 

Pd 



Tfr\' ■ (1) 



Here A {B + A) and a (b) denote the target (excited intermediate nucleus) and projectile (ejectile), respectively. The 
B -\- A system de-excites to the residual nucleus | (pc) by emission of the particles c and d which carry four- momenta 
{Ec,Pc) and {Ed,Pd), respectively, rm stands for the mass of particle i (i = A, C, a, 6, c, d) and (wc.e., Qce.) denotes the 
four-momentum transfer in the excitation process. In case that one of the outgoing particles is a boson (e.g. in the 
quasi-free decay) the according normalisation factor M/E has to be replaced by 1/2E. The full four body kinematics 
in the final reaction channel is included. 

The transition amplitude Tfi for the decay process is defined as 

Tft = {(fie, [MPc) <Pd{Pd)],^A) I ^cd,A I -0) (2) 

where Vcd,A denotes the A decay interaction that will be specified later. (j)c{Pc) and 4>d{pd) are the distorted wave 
functions of the outgoing particles c and d, respectively, and ipc is the wave function of the residual nucleus. The 
index [A) indicates that in case of the 2p emission the wave functions of the two outgoing identical fermions have to 
be antisymmetrized. The wave function | ij]) describes the intermediate B + A system and is defined by ||ll[ 

\^)^G\p)= \- I p) (3) 

bj + ii a/2 - Mb - J a - l/a - V\NAN 

where | p) is the doorway state excited initially by the reaction. The Green's function G describes the propagation 



of the {B + A) system and is approximated by that of the isobar-hole model [g^ 19,2C|,|2l|. rA(i^) is the energy 
dependent free decay width of the A, Hb is the Hamiltonian of the hole nucleus B, Ta and ?7a are the kinetic energy 
operator and the A-nucleus one-body potential, respectively, and VAAr,AAf is the A-hole residual interaction. For the 
calculation of we use the same input parameters as used in Refs. [p|Jll]. The Pauli blocking effects are assumed 
to be included in the average A-nucleus one-body potential, which we fixed by re-analysing the relevant scattering 
data H,^. We refer the reader for more details to Refs. [ppl]]. 

For charge exchange reactions the doorway state | p) has the explicit form 

I Pee.) = (xi^Vfc I tNN.NA \ xi^^Va^A) (4) 

where Xo^"* and xi denote the projectile and ejectile distorted wave functions, respectively, and (pa and tpb are the 
corresponding intrinsic wave functions of the projectile a and ejectile b; \ lpa) describes the target ground state wave 
function. The effective NN N A transition operator for the charge exchange process is denoted by tMN,NA- The 
round bra ( | on the right side of eq. (^ denotes the integration with respect to the projectile coordinates only. 
For pion induced reactions the coincidence cross section in the LAB frame is threefold differential 

d^cr _ 1 1 1 riic md mg PcEcPdEdEg ^ , ^ p 

dEcd^cdild ~ (2^)5 I Vrel \ 2uj^ E, Ed Ec I {mA E,) - ^{pa cosOad - Pc cosO^d) I ^ ' ' ^ ' 

Pd 

where (cj^r, ^-k) is the four- momentum of the mcident pion. Now the doorway state has the form Q 

^""^^ ^q^-S'^\ Tl e'«-^^|^A) (6) 



where S'f and ft are the spin and isospin transition operators, respectively, that convert a nucleon into a A(1232) 
isobar. The coupling constant f-nNA is fixed from pion-nucleon scattering data and has the value /^tva/^"" = 0.324. 
The index v — ±1 distinguishes between tt^ scattering. 
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FIG. 1. Schematic representation of the 2p emission induced by the decay of tlie A: (a) direct graph, (b) exchange graph. 
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B. The decay interactions 

In case of the "'^^C(7r+,7r+p)^^B and ^^C('^He,t7r+p)^^B reactions the decay interaction Vcd,A of eq. (0) is represented 



NA 



(7) 



Note that this interaction has no free parameter and is known from elastic pion scattering in the A resonance region. 
The expUcit formulas for the quasi-free decay of the A are given in appendix A. 

In case of the ^^C(7r"'",pp)^°B and ^^C(^He,tpp)^"B reactions the interaction for the process A + iV— > A^ + A^is 
described by a tt + p + 5' model Q : 



with 



VppM^',q') = V^iu;',q") + Vpiu;',q') + V5 

f-rrNNfirNA 



Vs = he 



9'na {^2 ■ Si^ (j2 ■ Tl^ 



(8) 



(9) 



In eq. lu' and q' are the energy and three- momentum transfer involved in the interaction A -I- A^ ^ A?^ -I- A", 
respectively. The interaction Vs is the so-called Landau-Migdal term. It describes the short range correlations 
for A -I- A^ N + N transitions. The special value for the Landau-Migdal parameter = 1/3 (in units of 

ficfTrNNf-KNA/rn^ ~ 800 MeV iir?) is known as the 'minimal (?5va' because it cancels out the attractive short range 
part of the 7r-echange potential [|3[ . 

The TT- and p-exchange potentials V.,^ and Vp are defined consistently with the potentials for the residual interaction: 



Vt,{uj' ,q) = tic 



firNNfirNA ( — ^ 



1 



le 



{a2 ■ q'){Si ■ q') T2 • f ] 



(10) 



Vp{uo\q') = nc 



fpNNfpNA 



A2 

p 



((72 X q') ■ {Si X q') (f2 • fi) 



(11) 



In the eqs. (|10| ) and (|l l|) , t' = lu —q is the four-momentum transfer in the decay process, and A^ (jUp, Ap) are the 
mass and cutoff mass of the tt (p), respectively. The various parameters are fixed as follows: fTiNNfTrNA/^TT = 0.162, 
fpNNfpNA/4:Tr = 8.32, = 0.14 GeV, mp =='^0.77 GeV, A^ = 1.20 GeV, and Ap = 2 GeV. 

As a consequence of the Pauli principle the wave functions of the two outgoing protons have to be antisymmetrized. 
This leads to two contributions to the 2p emission process, namely the direct and the exchange term (see Fig. 1). The 
antisymmetrization has to be carried out only for the finite-range tt- and p-exchange potentials. The Landau-Migdal 
term is a zero-range interaction and has not to be antisymmetrized. This treatment of the A -I- A^ ^ A^ -I- A^ interaction 
is in line with microscopic nuclear structure calculations |2^j2^ and microscopic G-Matrix calculations . Therefore 
the Landau-Migdal parameter g'j^/\ e xtracted from the 2p emission reactions can be directly compared with the values 
found in these calculations 23,2^]. The explicit formulas for the 2p emission process induced by the decay of the 
A in the nucleus are derived in appendix B. 
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FIG. 2. Triple differential cross section of the ^^C(7r+,p7r+)'^^B reaction atgiven pion and proton angle plotted versus the 
kinetic energy of the decay proton. The data have been taken from Ref. [E7|. The dashed curve shows the result without 
inclusion of the residual interaction Vajv,ajv while the dot-dashed curve shows the result with the inclusion of Vajv.ajv- The 
solid curve shows the result with additional inclusion of the distortion effect on the decay pion and proton wave functions. 



III. RESULTS AND DISCUSSION 



With the formalism described in Sec. II we have calculated cross sections for the quasi-free decay reactions 
i2C(7r+,7r+p)"B and i2c(3He,t7r+p)"B and the 2p emission reactions i2c(7r+,7r+pp)i°B and i2c(3He,t7r+pp)iOB. 
While the pion induced reactions are truely exclusive with the decay particles measured in coincidence at given energies 
and angles the decay cross sections of the charge exchange reaction have been integrated over a certain kinematical 
range, as determined by the geometry of the DIOGENE detector 10,1). Therefore the integrated cross sections for 
(■^He,t) induced processes are in a way less exclusive than the pion induced reactions and show only the gross features 
of the process. 

Furthermore, our calculations treat the distorsion effects on the incoming and outgoing particles in the adequate 
frameworks. In the pion induced reaction the distorsion of the incoming pion is treated within the isobar hole model 
while the distorsion of the decay pion and protons is described by optical model wave functions p5| , p6[ |. In the ('^He,t) 
charge exchange reaction the projectile and ejectile and the decay particles are described by optical model wave 
functions, whereas the A propagation through the nucleus is again treated within the isobar-hole model. 



A. The reaction ^^C(7r+,p7r+)"B 



In Fig. 2 we compare the results of our calculations for the ^^C(7r+,p7r+)^^B reaction to the data of Ref. p^ ] 
at the pion kinetic energy of T^r = 245 MeV. In the experiment the angles of the outgoing proton and pion were 
fixed at dp = 37.5° and 6*^ = —130°, respectively, in coplanar geometry. The threefold differential cross section is 
plotted versus the kinetic energy of the outgoing proton. Three different calculations are compared to the data. The 
dashed curve shows the result of our calculations with Van, an = while the dot-dashed curve shows the result with 
Van, AN 7^ 0. By comparison of both curves one recognizes that the inclusion of the residual interaction reduces the 
quasi-free decay cross section by a factor of ^ 4. This reduction is due to the absorption taking place in the multiple 
scattering of the pion. A similar reduction factor is also observed in the total pion- nucleus cross section The 
solid curve shows the result with additional inclusion of the distortion effect on the outgoing proton and pion wave 
functions. The distortion effect leads to a further reduction of the cross section by a factor of ~ 3. We describe the 
relative motion of the decay particles with respect to the residual nucleus by optical model wave functions. This is a 
consistent method within the framework of direct nuclear reaction theory and has been used in the analysis of other 
reactions, like A(e,e' pp). For the calculation of the proton and pion wave functions we used the optical potential 
parameters, as derived from elastic proton- nucleus [ p5[ and elastic pion- nucleus scattering p6[ |. Using these optical 
model wave functions we overestimate the absorption and thus obtain a lower limit for the~^^C(7r+,p7r+)^^B cross 
section. A comparison of the solid curve with the data shows that we underestimate the data by ~ 10 — 20% . This 
is in agreement with our expectation and assures us that our treatment of the distortion effects is reasonable. In the 
following calculations we will use the same model for the description of the distortion effects on the outgoing particles. 
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FIG. 3. Coincidence spectra for the ^■^C('^He,t) reaction at Tue = 2 GeV and triton scattering angle 9t =0°. (a) The 
'^^C(^He,tp7r+) data in comparison with the result of our calculation, (b) The ^^C(^H[e,tpp) reaction data The solid and 
dashed curve show the results with and without inclusion of the residual interaction, respectively. 




FIG. 4. Double differential cross section of the ^■^C(7r+,pp)^''B reaction. The angle of proton 1 is fixed at 6^1 = 130°. The 
cross section is plotted versus the angle of proton 2. The data have been taken from Ref. |28|p9|] . The dashed curve shows the 
result with Vajv,ajv = while the dot-dashed curve shows the result with Vaat.ajv 0. The solid curve shows the result with 
additional inclusion of the distortion of the decay protons. 



B. The reaction ^^C(7r+,pp)"'B 



In Fig. 3 we show the results for the "'^^C(7r"'",pp)^°B cross section at T^r = 165 MeV. The experiment was performed 
in coplanar geometry with the angle of proton 1 fixed at 9pi = 130°. The threefold differential cross section is plotted 
versus the angle of proton 2. The dashed curve shows the result of our calculation with Vaw.aat = while the dashed- 
dotted curve shows the result with VAAf.AW 7^ 0. The inclusion of the multiple scattering of the pion in the medium 
results again in a reduction of the cross section by a factor of ^ 4. The solid curve shows the result with additional 
inclusion of the distortion effect on the outgoing proton wave functions. This leads to a further reduction of the cross 
section by a factor of ~ 2. The calculations shown here were performed with a minimal Landau-Migdal parameter 
g'ffi^ = 1/3 in the 2p emission matrix element, which gives a very good agreement with the data |28|. Neglection of the 
TT- and p-exchange potentials in the decay interaction leads to almost the same results. Thus the Landau-Migdal term 
is the most important ingredient to the A-f — > N + N interaction. This is a consequence of the large momentum and 
energy transfer involved in the A + iV ^ N + N decay process leading to a very short ranged interaction. We remark 
that the determination of the gJvA parameter in the (tt^ ,pp)-reaction is rather direct since this reaction is dominated 
by the intermediate A++ excitation. On the other hand, in case of the ^^C(7r+,pn)'^°C reaction p3] or photon induced 
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2p emission reactions many competing processes are possible and the Landau-Migdal parameter g'l^^ can only 

be extracted in an indirect way. 

In Fig. 4 we study the same reaction as in Fig. 3; here the two protons were detected at Opi = 82.5° and 6p2 — —70°, 
and the cross section is plotted as function of the kinetic energy of proton 1. The curves shown have the same meaning 
as in Fig. 3. Note that the data are given in arbitrary units. Thus the data are normalized to the solid curve. The 
inclusion of the residual interaction and the inclusion of the distortion effects reduces the magnitude again by a factor 
of ~ 4 and ~ 2, respectively. The dependence of the cross section on the kinetic energy of the proton 1 is reproduced 
well. The slight deviation of the calculated cross section from the data at low kinetic energies could indicate that the 
treatment of the distortion on the outgoing nucleons should be improved for these energies. 



C. The reaction ^^C(^He,tp7r+) 

In Fig. 5 (a) we compare our results for the ^^C(^IIe,tp7r+) reaction with the data of Hennino et al. The data 
are integrated over the kinematically allowed pion and proton energies as well as pion and proton angles. The solid 
curve shows our calculation with inclusion of the residual interaction VAAr,AJV and with inclusion of the distortion 
effects on the outgoing particles. We see no effect of the residual interaction in this energy and angle integrated cross 
section, i. e. the cross sections with and without inclusion of Vaa^.a^v are the same. This is an indication that the 
p7r+ events in this reaction come only from the nuclear surface. 
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FIG. 5. Threefold differential cross section of the ^^C(7r^,pp)^''B reaction. The two protons were detected at O^i — 82.5° 
and 6p2 = —70°, respectively. The cross section is plotted as function of the kinetic energy of proton 1. The data have been 
taken from Ref. |2^,^. The curves shown have the same meaning as in Fig. 3. Note that the data are given in arbitrary units 
and are thus normalized to the solid curve. 



D. The reaction ^^Cf He,tpp) 

In Fig. 5 (b) we compare our microscopically calculated ^^C(^He,tpp) coincidence spectra with the data of Hennino 
et al. The data are integrated over the phase space of both outgoing protons. Two calculations are compared to the 
data: The solid curve shows the result with VAAr,AAf 7^ while the dashed curve shows the result with Vajva^^ = 
both cases the distortion effects are taken into account. In order to reproduce the magnitude of the experimental data 
we used a Landau-Migdal parameter of = 0.28. This value is consistent with that found in the (7r+,pp)-reactions 
of Figs. 3 and 4. Both values for g'f^^ lie in the range from 0.25 - 0.35 and are thus in agreement with values of the 
Landau-Migdal parameter found in microscopic G-Matrix calculations p3]. 
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IV. SUMMARY 



In conclusion, we have presented microscopic calculations for the quasi-free decay of the A resonance and the 2p 
emission in nuclei induced by pion absorption and by charge exchange reactions. These calculations are performed 
within the framework of the A-hole model and are consistent with our former calculations of inclusive and exclusive 
reactions on nuclei. Since the coupling interaction for the quasi-free decay is known we use this process to study the 
distortion effects on the outgoing pion and proton wave functions. In the 2p emission reactions we describe the A-f N— > 
N-|-N decay interaction by a n + p + g' model, which is consistent with our description of the residual interaction. 
We find that the A + N ^ N + N decay interaction is dominated by the zero-range Landau-Migdal term. The data 
for both the (7r+,pp) reaction and the (^He,tpp) reaction are well reproduced by calculations with a Landau-Migdal 
parameter in the range of (J^VA ~ 0.2b — 0.35. 

This work is supported in part by the U.S. Department of Energy under Contract DE-FG05-84-ER40145. 



APPENDIX A: EXPLICIT FORMULAS FOR THE QUASI-FREE DECAY 



In this appendix we derive the explicit formulas for the transition amplitude of the quasi-free decay. Inserting the 
decay interaction V^tt.a of eq. (Q) into eq. @) one can write the transition amplitude in the following way: 



(Al) 



Here Jj^ and represent the total angular momentum and the magnetic quantum number of the final (A-l)-nucleus; 
q^ and are the three-momentum and the spin projection of the outgoing proton; is the three-momentum of the 
outgoing pion. We expand the wave function | ?/') of the A-hole state in terms of the channel wave functions 

I [i^A^/iljtmJ = (jAmAjhmhljtmt) \Yj^rn^<S>j^^J , (A2) 

where YjA,mA the spin-angle wave function of the A and ^j^nih is the hole wave function of nucleus B. Thus we 
obtain 



Wo , 



(A3) 



jtmt Ah 



In ( [A^ ) Nc denotes the total number of allowed A-hole states. The radial wave function is then given by the inversion 
of eq. @: 



The wave functions of the decay nucleon and the outgoing pion are also expanded in multipoles: 



Ipmip jpfUp 



JplTlp 



(A4) 



(A5) 



(A6) 



where Xip ^^'^ Xi^r denote the distorted radial wave functions of the outgoing proton an d pion , re spectively. Using 
the expansion of eq. (A3) for the wave function | ?/;) and using the expansions of eqs. ( [A5| ) and ( |A6| ) for the outgoing 
nucleon and pion wave functions we can rewrite the transition amplitudes for the quasi-free decay of the A as 



T, 



q-f- 



jtmt Ah 



2V47r 



3 A I A ll 
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Yi^rm^ {%) (-l)-'""™" Yj^^^{q!^) (ja jh mi \ jt mt) (jA Jtt - | jp rUp) \ Ia ^ Ja 

[Ia 0) 1 I jr, 0) {Ip mi^ i m; | m^) (1 ^ i | f ta) J dr [q'^r) xi, Wpr) Xi. (q'^r) i^f^'^ (r) . 

(A7) 



In eq. (A7) the Clebsch-Gordan coefRcient (1/^ i | |ta) describes the isospin coupHng coefhcient of the decay 



process and a; — \j2x + 1. 



APPENDIX B: EXPLICIT FORMULAS FOR THE 2P EMISSION 

Here we show the exphcit formulas for the 2p emission processes. Due to the antisymmetrization of the two 
outgoing protons the transition amphtude for the 2p emission cons ists of the sum of the direct and exchange transition 
amphtudes. As a consequence of the arguments given in section II E we have to antisymmetrize only the tt- and p- 
exchange interactions. The transition amplitude is given by (see Fig. 1): 



1 



^2 V2 



pp,A 



x/2 



(414; 1; 2 I 14 + + I V-) - {J'kMk. 2-\\V.^Vp\ i,)) 



(Bl) 
(B2) 



In eqs. (Bl) and (B2) 4 ^^id -^^^ represent the total angular momentum and the magnetic quantum number of the 
final (A-2)-nucleus; q'^x and m^i ((7^2, "^52) are the three-momentum and the spin projection of the outgoing proton 
1 (2), respectively. Using the exp ansion of eq. (A3) for the wave function | -i/;) and the multipolc expansion of the 
nucleon wave functions of eq. ( |A5|) we can write the matrix element for the direct decay graph in the following way: 



(4M;;i;2|y,„Ai^)^EE YYY^^^(^-) 

jtmt Ah ipimi^i Jpi™pi 'p2™ip2 Ip'^'^p'^ ^h^^'h2 •■'p^^p •^i-'a 

,h+;2-;pi-v+;;2+^A (_i).A+./2-.. (_i)./2+Af2 jUiJpAj'h2yjp^jAhi2i^iA 

(jpi rupi jp2 mp2 I Jp Mp) {Jp Mp J'^ M/, | jt nit) Yi^i™,^^ (^pO ^p2m,^2 (9p0 
(/pi m;^! i m(,i I jpi mpi ) (/p2 m;^^ i m^2 | jp2 mp2) W{J2 Jpi jt jh ; JA Ji) 

Jp^ J2 

Jp J'h jt . 
1 

("^) E (1 5 V I 1 ^a) (1 m i ^2 I i Tp2) 

m— — 1 



/pi 


1 

2 






/p2 


1 

2 


Jp2 ^ 


Za 


3 
2 


JA 


II 


^'h2 


1 
2 


fh2 


^1 


1 


J2 . 




. I2 


1 





Here 



dr2 rl <pn'^^ (V^^ i) j;, {^2) X/p2 ('?p2''2) j dn n vi,^i^j^{ri,r2) xi^i iq'p^ri) V'(a/o (^1) • 

(B3) 

denotes the radial hole wave function with quantum numbers nj^a /^2 i^2 . The non-local potential 



'!^;i,Z2,J2(''i;^2) is the sum of the central and tensor part of the decay interaction |35 
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VhMAri,r2) = -^ir 5i,^jji,.j, Vgir^.r^) + (-if- hh {h(^l2Q\ 20) W^i 1^2 1 ; J2 2) ^z,,,/, (^i, ra) . 

(B4) 

For the direct matrix element Vj^ and V^^^ are the multipole expanded central and tensor parts of + + Vs [|2| ; 
for the exchange matrix element they are the multipole expanded central and tensor parts of V-^ + Vp. 
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